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Abstract. A systematic study of the production of cj-mesons in proton-proton-collisions was carried out 
in a kinematically complete experiment at three excess energies {e = 92, 128, 173 MeV). Both protons 
were detected using the large-acceptance COSY-TOF spectrometer at an external beam line at the Cooler 
Synchrotron COSY at Forschungszentrum Jiilich. The total cross section, angular distributions of both 
a;-mesons and protons were measured and presented in various reference frames such as the overall CMS, 
helicity and Jackson frame. In addition, the orientation of the tj-spin and invariant mass spectra were 
determined. We observe w-production to take place dominantly in Ss and Sp final states at e = 92, 128 
MeV and, additionally, in Sd at e = 173 MeV. No obvious indication of resonant a;-production via A'^*- 
resonances was found, as proton angular distributions are almost isotropic and invariant mass spectra 
are compatible with phase space distributions. A dominant role of ^Pi and ^So initial partial waves for 
w-production was concluded from the orientation of the decay plane of the tj-meson. Although the Jackson 
angle distributions in the ojp- Jackson frame are anisotropic we argue that this is not an indication of 
a resonance but rather a kinematical effect reflecting the anisotropy of the lu angular distribution. The 
helicity angle distribution in the a;p-hehcity frame shows an anisotropy which probably reflects effects of 
the CO angular momenta in the final state; this observable may be, in addition to the orientation of the 
w decay plane, the most sensitive one to judge the validity of theoretical descriptions of the production 
process. 

PACS. 13.75Cs Nucleon-nucleon interactions - 13.88e-)- Polarization in interactions and scattering - 
14.40Cs Other mesons with S=C=0, mass < 2.5 GeV - 25.10+s Nuclear reactions involving few-nucleon 
systems - 25.40Ve Other reactions above meson production thresholds (energies > 400 MeV) 



1 Introduction 

The production of w-mesons in proton-proton collisions 
has come into the focus of experimental and theoretical 
studies only in recent years. Early measurements of pro- 
duction cross sections employed hydrogen bubble cham- 
bers, a compilation of those data was published by Flaminio 

Correspondence to: w.ullrich@physik.tu-dresden.de (W. Ull- 
rich) 



et al. 1 for beam momenta of pbcam = 3.99 to 19.0 GeV/c. 
In the last 10 years several experiments have been per- 
formed with electronic detectors, namely by the DISTO 
collaboration at Pbcam = 3.67 GeV/c [2], [3] and at beam 
momenta below 3.2 GeV/c by Hibou et al. [3], Barsov et 
al. [5], and the COSY-TOF collaboration [5], [7|. These 
experimental data provide, in principle, important infor- 
mation about the short range part of the nucleon-nucleon 
interaction which is dominated by the exchange of the 
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Fig. 1. Tree level diagrams for the reaction pp — >■ ppcj. Emis- 
sion from (a) nucleon line, (b) internal meson line, and (c) 7V*- 
resonance. Initial and final state interaction is not indicated. 



isoscalar oj-meson They also furnish elementary cross 
sections needed for studies of in-medium properties of 
mesons. 

Various theoretical models were developed to describe 
the experimental data. Sibirtsev [3] has investigated the 
p,uj, and (/)-production within a one-pion exchange model. 
Nakayama et al. |10j described the w-production within a 
meson exchange model taking into account nucleonic and 
meson-exchange currents. These studies were extended in 
[llj and [12. . The w- production via nucleon resonances 
was studied by Fuchs et al. [TB] and Faessler et al. [TJ. 
Studies which include both resonant production and me- 
son currents were performed by Kampfer et al. [15j and 
Titov et al. [16j . The diagrams considered so far in theo- 
retical calculations are shown in fig. [TJ 

These models were particularly tailored to reproduce 
the total cross sections and/or angular distributions of the 
w-meson. The rather limited experimental data basis does 
not suffice to prove one or the other model preferable. In 
addition, none of the theoretical models takes into account 
all diagrams of fig.[l]and most of them neglect either initial 
or final state interactions or both. 

In this paper we present a wealth of experimental data 
for the pp — >■ ppw reaction obtained at three different 
beam momenta, namely pbeam = 2.950, 3.059, 3.200 GeV/c, 
which correspond to excess energies (e — y/s — 2mp — m^) 
of 92, 128, and 173 MeV, respectively. All data were taken 
with a special downscaled trigger in parallel to studies of 




Fig. 2. Sketch of the COSY-TOF spectrometer with its main 
components; see text for details. 

associated strangeness procuction. First evaluations of the 
data at 2.950 and 3.200 GeV/c, taken in 2000, focused 
on the angular distributions and total cross section for 
w-mesons 0. Here we show data where the method of 
analysis was improved and extended to other observables. 
The data at 3.059 GeV/c were collected in an experiment 
carried out with a very large integrated luminosity in 2004 
as it was devoted to the search for a supposed pentaquark 
state jTT]. Other publications of our collaboration were 
concerned with the measurement of the analysing power in 
the same reaction channel using a polarised proton beam 
at 3.065 GeV/c [TH] and a preferred orientation of the w 
decay plane [19]. 

This paper is organised as follows: The experimental 
procedure is detailed in sect. 2, describing the detector 
setup, the principle of measurement and data analysis, 
Monte Carlo simulations, and a discussion of systematic 
uncertainties. In sect. 3 we present and discuss the ex- 
perimental results obtained for the total cross sections, 
angular distributions of the w-mesons and protons in the 
overall centre-of-mass system (CMS), the orientation of 
the Lo decay plane, invariant mass spectra, and distribu- 
tions of both helicity and Jackson angle in their respective 
frames. These results set important benchmarks for any 
theoretical model. The paper ends with a summary. 



2 Experimental procedure 

2.1 Detector setup 

The experiments were carried out with the time-of-flight 
spectrometer COSY-TOF installed at an external beam 
line of the COoler SYnchrotron COSY at Forschungszen- 
trum Jiilich. A sketch of the detector is shown in fig. [21 

The extracted proton beam (spill length « 5 min, sev- 
eral 10^ protons/s) is directed onto a liquid hydrogen tar- 
get of 6 mm diameter and 4 mm length [20] through an- 
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Fig. 3. Monte Carlo data for pp— )-ppa; (cj 



"'"tt'tt'^), pp— >-pp7r'''7r~7r'' and pp— >pp7r"''7r~ at e = 128 MeV. The experimental 
resolution of the detector has been accounted for. Three regions (A, B, C) for the separation of protons and pions are indicated. 
See text for details. 



nular veto detectors of different diameter positioned at 
260, 51 and 3.5 cm upstream of the target. The emerging 
particles traverse a layered system of time-of-flight start 
and tracking detectors (called Erlangen start [21]). After 
a flight path of ~ 3 m through the evacuated tank (0.2 
Pa) all charged particles are detected in the highly granu- 
lar stop components of the spectrometer. They consist of 
two three-layered forward hodoscopes (quirl, ring) which 
are built according to the layout described in [32, and 
a barrel hodoscope [53], all utilising BC412 scintillating 
material. 

From the combined measurement of time and position 
the velocity vectors of all charged particles are determined 
with a time-of-flight resolution of better than cttof — 300 
ps and an angular track-resolution of better than — 
0.3°. Due to the low mass area density of the time-of- 
flight start and tracking detectors, the influence of small 
angle scattering and energy loss is almost negligible for 
particles with /3 > 0.5. Only particles with these velocities 
are produced in the reaction under study. 

Unlike magnetic spectrometers, which often provide 
particle identiflcation at the cost of limited acceptance, 
the COSY-TOF spectrometer stands out for its large geo- 
metrical acceptance (1° < 9ii,b < 60°, 0° < < 360°) and 
an efficiency > 95% for the detection of a charged parti- 
cle. This allows the almost unambiguous and simultaneous 
identiflcation of different reaction channels (e.g., pp— > pp, 
d7r+, ppTT+TT", ppw, pK+A, pK+I]°, pK"S+) by exploit- 
ing the time of flight of all emerging charged particles as 
well as the event's topology. 



2.2 Principle of measurement and data analysis 

The reaction pp — > ppw with the main decay channel of 
the w-meson (w tt+ tt- 7t° , BTZ « 89.2% [^) is prese- 
lected during data taking via a trigger set on four charged 
particles detected in the stop components (quirl, ring, bar- 
rel). The outgoing protons and the charged pions are, due 
to their large mass difference, clearly separated in a veloc- 



ity versus polar angle representation as can be seen from 
the left frame of fig. [31 cf. also [7]. 

This kinematical separation serves to identify protons 
as they are restricted to a region with 9 < 23° and mod- 
erate velocities, while the pions mainly cluster at /3 w 1 
over the full angular range. There is only a small number 
of events (about 10 to 15% depending on excess energy) 
where pions are found in the proton region. For the anal- 
ysis we require exactly two particles within one of the 
indicated proton regions A, B, or C, shown in the left 
frame of fig. [3] (the effect of these different regions will 
be discussed below), and exactly two particles outside, 
which are, in turn, assumed to be pions. In addition, all 
tracks must fall into the region of optimum geometrical 
acceptance (3° < 6* < 60°). Proton and pion masses are 
then assigned accordingly and the four-momenta of pro- 
tons are calculated using the measured velocity vectors. 
Monte Carlo simulations showed that the assignment of 
protons and pions to the pp— J'ppo; reaction channel based 
on these criteria is correct for 99.2% of all events. We 
would like to point out that the geometrical acceptance 
for pp — >■ ppoj — >■ pp7r"'"7r~7r'^ reactions, that means the 
fraction of events with all 4 charged particles in the geo- 
metrical covered region, is 65% at the lower excess energy 
£ = 92MeV and increases to 68% at e = 173MeV. This 
applies for phase space and isotropic distributed events. 

The missing mass distribution calculated from the mo- 
menta of the identified protons exhibits a peak at the w- 
mass. However, this spectrum contains a large contribu- 
tion of resonant two pion production via pp— >■ ppp, p — > 
TT+TT^ and non-resonant production of two pions (pp — >■ 
pp7r+7r~) as well as three pions (pp-^pp7r+7r~7r'^), all of 
them constituting the major components of an unavoid- 
able physical background. In order to accentuate this point, 
we show in the middle (right hand) frame of fig. [3] Monte- 
Carlo simulated data for non-resonant 7r"'"7r~7r° {tt~^'it~) 
production. It is obvious that in these two cases the kine- 
matical separation of protons and pions is less certain. 

It is possible to reduce the two-pion background by 
exploiting the topology of the w-decay into three pions. 
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Fig. 4. Proton-proton missing mass spectra at e = 128 MeV 
obtained when applying the different regions for proton assign- 
ment shown in fig. |31 



Here, the plane defined by the two charged pions does, 
in general, not contain the pp-missing momentum vec- 
tor (= w-momcntum vector) due to the momentum of 
the undetected tt*^. By investigating Monte Carlo data 
it was found that requiring an acoplanarity angle a — 
ACfn^ X -f^,),{t., X 1^^)) > 5° rejects 90% of the 
two-pion background while only 17% of the w-events are 
removed, thus improving the signal-to-background ratio. 
We would like to stress the robustness of the method as the 
selection of different acoplanarity angles result in the same 
yield of the w-signal after acceptance correction, within an 
uncertainty of 2% [7 . Throughout the analysis an acopla- 
narity angle of a > 5° is required. Finally, in the anal- 
ysis only those events are taken into account for which 
the sum of the proton momenta points into the backward 
hemisphere of the overall CMS in order to improve on 
the momentum resolution. Because of the identical parti- 
cles in the entrance channel all distributions of final state 
particles in the overall CMS will be symmetric so that a 
reduction of information does not ensue. 

In the present paper our former study [7] is buttressed 
by scrutinising the influence of the proton region chosen 
on the yield determined for the w-signal. As already men- 
tioned, three proton regions are considered, region A with 
rather narrow boundaries, a wider region B (formerly used 
in [7]) and an even larger one, region C. As shown in 
fig. m a distinctive peak at the w-pole mass (m^^ =782.7 
MeV/c^, [21] ) can be seen in all cases above the multi- 
pion background, the shape of which for missing masses 
smaller than the w-mass depends dramatically on the re- 
gion chosen for the separation of protons and pions, while 
the shape above the w-mass is barely different. 

Region A, which sharply restricts the protons in angle 
and velocity, results in a missing mass spectrum where the 
shape of the background can hardly be estimated in the 
region of the w-mass. The missing mass spectrum result- 
ing from region B exhibits an w-signal which resides on 
a smooth, apparently convex background which extends 



to missing masses smaller than the oj-mass. When further 
enlarging the proton region to region C the background 
dominates the spectrum and its shape below the cj-signal 
tends to change from a convex to an almost linear shape. 
The advantage of using region B [7] lies in the best signal- 
to-background ratio, while in case of region C the w-signal 
becomes less susceptible to the method used for its extrac- 
tion (see below). This results in a smaller uncertainty for 
the extracted counts, although about 20% of the signal 
is lost. Thus, in contrast to our earlier analyses [7] we 
now use the wide region C which is adjusted with increas- 
ing excess energy in order to accommodate the protons 
which become less confined in angle and velocity. It is in- 
teresting to note that a signal of the 77-meson appears at 
550 MeV/c^, however, the excess energy of about 350 MeV 
in this channel causes a large width of this signal greatly 
complicating a detailed analysis of the data. 

In order to quantitatively determine the yield of the 
w-signal above background we follow two different ap- 
proaches, briefly described in [19j. One method is a simul- 
taneous least square fitting of a Voigt function (convolu- 
tion of a Gaussian and a Breit-Wigner function) for the w- 
signal and a second-order polynomial for the background 
below the signal. The width of the Breit-Wigner distribu- 
tion is fixed to the natural width of the w-meson (F = 8.49 
MeV, [Ml), while all other parameters are allowed to vary 
freely. The alternative method is based on the simulation 
of the main background channels (pp — ppp, p — tt+tt", 
pp ppvr+TT^, pp — )• ppTr+vr^vr*') where the simulated 
data are analysed in the very same way as the experimen- 
tal data. The yield of various background contributions 
obtained are adjusted in order to reproduce the experi- 
mental background. It is found that the 37r-background 
dominates over the 27r-background (which was effectively 
reduced by the acoplanarity requirement); the non reso- 
nant 27r-contribution can be exchanged for the resonant 
one from p-decay as the p-meson has a large width. The 
result of both procedures is shown in fig.[S]for an excess en- 
ergy of 128 MeV. The background described by the "pion- 
cocktail" agrees very well with the one determined by the 
fitting procedure. The yield determined for the w-signal 
is found to agree within a few percent for the two meth- 
ods; the present results are fully consistent with but more 
precise than our earlier ones. Although this experiment 
does not aim at a measurement of the w mass we like to 
quote the parameters determined for the Voigt function at 
the three excess energies of £ = 92, 128, 173 MeV, namely 
= 780.2, 782.4, 779.5 MeV/c^ with a systematic un- 
certainty of 0.3% and a standard deviation of the Gaus- 
sian resolution function of a = 7.6 ± 1.1, 10.3 ± 0.5, 11.3 ± 
0.6MeV/c^. As usual, the best momentum and, hence, 
mass resolution for time-of-flight experiments is achieved 
for the slowest particles, i.e. protons at the smallest excess 
energy. 

The procedure of yield determination must be carried 
through for each interval of the differential observables, 
e.g. for different intervals in cos 9^ to obtain the u! angular 
distribution in the overall CMS. 
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Fig. 5. Proton-proton missing mass spectra at e = 128 MeV 
with adjusted Monte Carlo distributions obtained for resonant 
pp— >-ppp, p tt'^-k' production and non-resonant two- and 
three-pion production. The sum describes the background be- 
low the oj-signal very well. (See text for details.) 



The absolute normalisation is accomplished by evalu- 
ating elastic scattering, which is always measured simul- 
taneously at COSY-TOF. The angular distribution deter- 
mined is compared to literature data j25| . whereupon the 
normalisation factor directly yields the integrated lumi- 
nosity. The uncertainty of this procedure (5%) is in equal 
parts due to the intrinsic uncertainty of our measurement 
and the error of the literature data. 



2.3 Monte Carlo simulations 

During the analysis of the experimental data each step 
of the analysis was checked against results obtained by 
Monte Carlo simulations, thereby giving insight into the 
detector performance and, in particular, its acceptance. 
The programme package [26','27' models the whole COSY- 
TOF-detector in great detail (geometry of all sub-detectors, 
detector housings and wrappings, cabling, supporting struc- 
tures etc.). The particles are generated according to equal 
population of phase space |28] and propagated in space 
and time accounting for small angle scattering, nuclear re- 
actions, 5-electrons and energy loss in all detector compo- 
nents; the energy loss value is used to generate a charge or 
light output signal. The AD-conversion takes into account 
light attenuation in the hodoscopes, efficiency and gain of 
the individual photomultipliers. The time and charge reso- 
lution of the individual elements of different sub-detectors 
is also accounted for. A properly chosen amount of noise 
is added to each signal In addition, various decay 
anisotropies, e.g. that of the w-decay (J^ = 1~) can be 
included optionally [50] • The Monte Carlo simulations also 
account for the anisotropic cos 0* angular distributions as 



experimentally observed (see section 3.2.1) at the three 
excess energies. These Monte-Carlo data are subjected to 
the very same analysis routines as the experimental data 
in order to treat both on an equal footing, thus, for all 
observables the detector acceptance can be determined. It 
is found that the acceptance is a smooth function for all 
differential observables as can be inferred from the figs. 
[7] to [15] shown below. The numerical value of the integral 
acceptance, however, depends largely on the size of re- 
gion C chosen for proton identification and turns out to 
be slightly different for the three excess energies studied 
(see table [T]). 



2.4 Systematic uncertainties 

Three major sources of systematic errors were identified 
for the present analysis. An obvious one is the uncertainty 
of the luminosity determination of 5% (see above). The 
second source of systematic uncertainties stems from the 
determination of the yield of the w-signal residing on a 
sizable background. The acceptance correction itself is the 
third source of systematic uncertainties. The Monte Carlo 
simulation should model, besides other observables, the 
relation /? vs 9\ah as perfectly as possible since it is used 
for particle identification. Although the majority of pious 
cluster at /3 ~ 1 some happen to fall into the /3 vs 6'iab- 
region chosen for proton assignment. This fraction is fairly 
sensitive to the correct simulation of the time resolution 
for pion tracks. Its infiuence can be studied by varying the 
/3-boundary of the proton region. Thereby the systematic 
uncertainty of the acceptance correction was determined 
to amount to 4% at e = 92 MeV and e = 128 MeV and 
6% at £ = 173 MeV. 

In the case of the total cross section these three sys- 
tematic uncertainties were added quadratically. 

In case of differential cross sections the systematic un- 
certainties for the determination of the yield of the uj- 
signal and the acceptance correction were investigated in 
detail for each bin of the observable under consideration. 
A spread of uncertainties was quantified, firstly, by a sys- 
tematic variation of the range of the missing mass included 
when fitting the background below and above the w-signal 
and, secondly, by a comparison of the results obtained by 
the two different fitting methods. The spread was found 
to be between 5 and 20% and was added to the statis- 
tical error. The sum is indicated as a common error bar 
for each interval of the observable under consideration. 
The systematic error due to the luminosity determination 
contributes an overall uncertainty which is not shown. 



3 Discussion of results 
3.1 Total cross section 

The number of w-events extracted from the proton-proton 
missing mass spectrum at different excess energies are 
used to calculate the total cross sections listed in table [T] 
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Table 1. Total cross section for the reaction pp — > ppoj at 
three excess energies; A'^^: number of cj-events; acc: detector 
acceptance as determined from Monte Carlo simulations; see 
text for details. 

e (MeV) iV„ acc / Ldt crtot (fJ-b) 

92 2052 0.156 1449 9.1 ± 0.6 ± 1.0 

128 116315 0.168 53500 12.9 ± 0.2 ± 1.1 

173 3558 0.155 795 28.9 ± 1.9 ± 3.1 



Table 2. Coefficients of Legendre polynomials determined by 
fitting eg.® to the u angular distributions. 



e (MeV) 




a'2 




92 


9.3 ± 1.1 


0.23 ±0.26 




128 


12.7 ±0.5 


0.49 ±0.10 




173 


30.5 ±1.8 


0.46 ±0.15 


0.42 ±0.18 



pointing into the backward hemisphere, hence we can show 
only the angular distribution for the forward hemisphere. 
A transition is observed from a slightly anisotropic to a 
pronounced anisotropic angular distribution. 

The standard parametrisation of an angular distribu- 
tion in the overall CMS in terms of Legendre polynomials 
reduces to 

^^Y. «2L-P2L, i = 0,1,2,... (1) 

since only even Legendre polynomials P2L need to be taken 
into account in the case of identical particles in the en- 
trance channel. This equation can be rearranged to 

^-^^d+Eoi.^/'..), (2) 

L=l 

where crtot /47r = qq and = a2L/o,o- We prefer the use 
of eq.lO over eq.([T|) as the former yields the total cross 
section crtot without further calculation being necessary 
and the significance of the contributions of higher partial 
waves in the exit channel is better judged via the ratio 
rather than a2L- The coefficients obtained are listed 
in tabled 

While it is sufficient to use Lj^nx — 1 for the two lower 
excess energies, as can be deduced from table [5] the inclu- 
sion of -Lmax = 2 improves the description of the angular 
distribution for e=173 MeV although the uncertainty of 04 
is rather large. The values of crtot are consistent, within the 
uncertainties, with those listed in table [TJ We like to point 
out that the deduced values of Lmax are fully consistent 
with kinematic considerations: the maximal orbital angu- 
lar momentum for the (pp)-a; system can be estimated as 
V = P* /"tLtt where p* is the centre of mass momentum in 
the final state, 1/(2to^) the strong interaction radius and 
TO^ the pion mass |14j . For the three excess energies rj 
varies from 2.4 to 3.3 h. However, the angular momentum 
barrier truncates this maximum orbital angular momen- 
tum to an effective one of L = l?i for e = 92 and 128 MeV, 
while L = 2h is only probable at e = 173 MeV. 

These results can be looked at from a different point 
of view. Meson production in pp reactions is usually de- 
scribed by defining the non-relativistic momenta p (rela- 
tive momentum of the two protons) and q (momentum of 
the w-meson with respect to the centre of mass of the two 
protons) and the associated orbital angular momentum 
quantum numbers L and l^i, respectively. Then the spec- 
troscopic notation (^'^+^Lj)i holds for the partial waves in 
the initial state and (^^~^^Lj)f l^^ for the final state which 
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Fig. 6. Total cross section as a function of excess energy mea- 
sured for the reaction pp — >■ ppoj. 



The values at e = 92 MeV and 173 MeV agree with those 
published earlier |5] , [7] , however with improved accuracy. 

Fig. E] shows all published data for this reaction; the 
present data follow the general trend. The results of the- 
oretical models more or less follow the data, but the to- 
tal cross section is not a sensitive observable that proves 
one or the other model preferable - additional observables 
are needed. As COSY-TOF almost completely covers the 
phase space of the reaction many differential cross sec- 
tion can be determined offering a variety of other physi- 
cal quantities setting benchmarks for theoretical models. 
Hence we turn to differential distributions. 



3.2 Angular distributions of the a;-mesons and protons 
in the overall CMS 

3.2.1 cij-mesons 



The uj angular distributions in the overall CMS are shown 
in fig. [7] for the three excess energies, the acceptance is 
shown in the lower part of each frame. It is a smooth 
function of cos 0^ which drops with increasing angle only 
slightly from «40 to «30%. As pointed out above, the data 
analysis requires two protons with an added momentum 
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Fig. 7. Angular distribution of the cj-meson in the overall CMS (upper part in each frame). The lower part shows the acceptance 
of the COSY-TOF spectrometer as determined from Monte Carlo simulations. 



Table 3. List of all partial waves in the initial state which can 
produce final states with two units of angular momentum at 
the most. The notation is defined in the text. 



type 


^2S + l 




amplitude 


Ss 




'Sos 


h 


Sp 


'So- 
'D2 - 


^ 'Sop 
'Sop 


h 

h 


Ps 


'So- 
'D2 - 
'D2 - 


>^PlS 
^PlS 
-+^P2S 


h 
h 
h 


Sd 


^P2- 
3F3_ 


^ 'Sod 
-> 'Sod 
'Sod 


fr 
fs 

h 


Pp 


^Po- 

^Pl^ 
^P2- 

spa - 

3F3_ 
^F4- 


^Po,l,2P 
^ ^Po,l,2P 
^Po,l,2P 
^Po,l,2P 
■>'Pl,2P 
■+'P2P 




Ds 


^P2- 

^Fa - 

^F3- 


-> 'Das 
'Das 
^ 'Das 


/25 

ha 
,(27 



in short will be designated as type "LIJ^ later on; S, L, J 
refer to the quantum numbers of spin, orbital angular mo- 
mentum, and total angular momentum of the pp-system 
in either initial (i) or final (f) state (see e.g. Hanhart [5T|). 
We list in table [3] all partial waves of the entrance chan- 
nel (antisymmetric pp-wave functions) which can produce 
(accounting for conservation of both parity and angular 
momentum) final states with at most two units of orbital 
angular momentum as suggested by the angular distribu- 
tion of the w-meson. 

A comparison with the value of Lmax deduced from 
the angular distribution shows that at £=92 and 128 MeV 
final states of type Ss, Sp and Ps (amplitudes fi,f2, --ife) 
contribute, while at e=173 MeV the final states of type 
Sd, Pp, and Ds (amplitudes f7,...,f27) come into play. It 
is important to realise that from a particular initial state 
two or more final states can be populated (cf. amplitudes 



Table 4. Coefficients of Legendre polynomials determined by 
fitting eq. ([2]) to the proton angular distributions. 



e (MeV) 




/ 

0.2 


92 


8.5 ±0.7 


0.11 ±0.22 


128 


12.7 ±0.3 


0.12 ±0.06 


173 


27.8 ± 1.3 


0.13 ±0.11 



/2 and /4 or /a, /s, and /g). Which of these amplitudes are 
present in the reaction and contribute with which weight 
depends on the transition matrix elements; they can be 
determined by a full model calculation. However, we con- 
jecture that, due to the available excess energy, the two 
final state protons are most probably in the ^Sq state. 
This would reduce the mainly contributing final states to 
those of type Ss, Sp, and Sd. It cannot be excluded that 
final states of type Ps are present as they originate from 
the same initial state as Sp type states, but probably with 
only little weight. 

3.2.2 Protons 

In fig. [5] we show the angular distributions of protons in 
the overall CMS for the three excess energies. The ac- 
ceptance (lower part of each frame) is a smooth function 
of cos 0* with an almost constant slope. The decrease to- 
wards cos 0* = 1 is due to the analysis requirement of the 
sum of the proton momenta pointing into the backward 
direction. The angular distributions are, within uncertain- 
ties, symmetric with respect to cos 9* = 0, as they should 
be in the case of identical particles in the entrance channel. 
This is a strong indication that the acceptance corrections 
were performed properly. 

Using even Legendre polynomials (Pq and P2) accord- 
ing to eq. the least square fitting resulted in the values 
of table m The distributions are almost isotropic; little 
room is left for an i = 1 partial wave. As a consequence, 
the final state pp-system is most probably in the ^Sg-state 
and less probably in the '^P-state; this conclusion is con- 
sistent with that drawn from the uj angular distribution. 

As mentioned in the introduction, there are theoreti- 
cal models which favour w-production to proceed via A^*- 
resonances (pp— >-pA^*, N* — pw). The finding of an 
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Fig. 8. Angular distribution of the protons in the overall CMS (upper part in each frame). The lower part depicts the acceptance 
of the COSY-TOF spectrometer as determined from Monte Carlo simulations. 



isotropic, at most slightly anisotropic proton angular dis- 
tributions at all energies and that of undeniably anisotropic 
ones for w-mesons at the two higher excess energies are dif- 
ficult to reconcile with resonant cj-production as the major 
reaction mechanism. This will be shown in the following 
discussion. 

In the present experiment the mass region of 1.721 
GeV/c^ (threshold) to 1.893 GeV/c^ (highest beam mo- 
mentum) was covered. In this region nucleon resonances 
are known which have significant N* pw decay branches, 
namely N(1710)Pii, N(1720)Pi3, and N(1900)Pi3 UM- Be- 
low the threshold N(1535)Sii, N(1650)Sii, N(1675)Di5, 
N(1680)Fi5, and N(1700)Di3 are listed in [24]. Most of 
them are classified as four star resonances, while N(1700) 
and N(1710) have three and N(1900) two stars. Above the 
accessible mass range the two star resonances N(1990)Fi7, 
N(2000)Fi5, N(2080)Di3 are located. All these N* res- 
onances have large widths of typically 100 - 200 MeV 
and, hence, their Breit-Wigner tail may well extend into 
the mass interval covered by this experiment. The in- 
clusion of subthreshold resonances has proven important 
in the calculation of coupling constants of baryon reso- 
nances to the nucleon-omega channel in particular the sub- 
threshold resonances Sn and D13; above the cj-threshold 
N(2080)Di3 was found to be of particular importance [35] • 
It should also be stressed that the N(1535)Sii resonance 
was specifically claimed to play a crucial role in off-shell 
w-production [T3]. Hence, in order to be as unbiased as 
possible, we include all of these resonances in the follow- 
ing consideration on resonant w-production. 

In the case of the pp —J- pN* , N* — )■ pw reaction the 
CMS angular distribution of the proton associated with 
the A^*-resonance depends on spin and parity of the A^*- 
resonance and can be calculated using the formalism de- 
veloped for nuclear reactions by Blatt and Biedenharn 
|33j . They showed that the transition from an initial state 
with J = Li + Si to a final state with J = Lf + Sf 
(here the final state is the pA^* system) is characterized 
by an angular distribution which can be described by an 
expansion in even Legendre polynomials Pl in analogy 
to eq. (H]). The expansion coefficients are uniquely deter- 
mined for a given set of quantum numbers Li, Si, J, L f, S f 
from Wigner and Clebsch-Gordan coefficients. The maxi- 
mum of the summation index is Lmax which, in the case 



of reactions with integer values of J, is given by £max= 
min(Li, J, Lf). Examples of angular distributions for pro- 
tons from the pp — ^ pA^*(1720) reaction for various sets of 
quantum numbers {Li, Si, J, Lf, Sf) can be found in |14| . 

We used this formalism 33 to find those pp — > pA^*- 
reactions which show isotropic proton angular distribu- 
tions. They are obtained only for transitions from the 
initial states ^Sq or ^Fq to any A^*-resonance and from 
■^Pi to Sii. This is the consequence of either Li, J 01 Lf 
being zero, thus Lmax = 0. Then, whatever resonance is 
involved, the angular distribution of the decay products is 
isotropic. Thus we have pinned down three initial states 
which can lead via resonant w-production to both isotropic 
proton and uj angular distributions. 

All other initial states result in angular distributions 
of various degree of anisotropy in the production step. In 
particular the transitions from initial states '^Pi and 
(amplitudes /i and /y of table |3]) to any A^*-resonance 
(except ^Pi to Sii which is isotropic) give angular distri- 
butions of the type oq + a2P2- If the iV*-resonance now 
decays isotropically, then the anisotropy from the produc- 
tion process is diluted. As A^*-resonances such as Pn and 
P13 decay with a relative angular momentum of Ifi, the 
anisotropy from the production process may be enhanced. 
This applies to the resonances N(1710)Pii, N(1720)Pi3, 
and N(1900)Pi3 in the mass region accessible in this ex- 
periment. We conclude that resonant w-production can 
produce isotropic as well as anisotropic angular distribu- 
tions for protons and cj-mesons, however the degree of 
anisotropy at e = 92MeV indicates that this contribution 
is of minor weight. It is also conceivable that resonant 
production provides a minor contribution at the higher 
energies. 

The initial states not mentioned explicitly so far yield 
proton angular distributions from the production step that 
are of the type qq + 02^2 + a4P4 and may even include an 
oePe term. This implies also an anisotropic sequential A^*- 
decay, and hence, anisotropic w angular distributions are 
expected. We are not aware of any mechanism that would 
yield isotropic proton angular distributions when superim- 
posing those of anisotropic production and anisotropic se- 
quential decay. Hence we conclude that the isotropic pro- 
ton angular distributions and the strongly anisotropic w 
angular distributions observed at the higher excess ener- 
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gies are incompatible with resonant tj-production domi- 
nating the reaction mechanism. 

We summarise the conclusions drawn from the inter- 
pretations of the angular distributions of cj-mesons and 
protons in the overall CMS: the former are produced with 
exit channel angular momenta of 0, 1, and 2 h depending 
on excess energy while protons are found to be emitted al- 
most isotropically at all excess energies. This implies that 
protons are produced most probably in the ^So (ampli- 
tudes fi, f2, fs, fy, fs, fg of table [31), and less probably in the 

final state (amplitudes f^, /s, /g). The contribution of 
the other amplitudes is assumed to be insignificant. The 
definite anisotropy of the w angular distribution and the 
almost isotropic proton angular distribution observed at 
e = 128 and 173 MeV rules out resonant w-production 
to be the major reaction mechanism at these excess en- 
ergies. It cannot be excluded that it contributes to the 
angular distributions of oj-mesons and protons to a minor 
extent at all excess energies as resonant cj-production via 
the A^*-resonances Sn, D13 yields isotropic and via Pn or 
P13 anisotropic angular distributions. In these cases, only 
particular initial states contribute. 



3.3 Orientation of the w-spin 

Already in 1962 it has been pointed out by Gottfried 
and Jackson |34l that the orientation of the w-spin gives 
valuable information on the reaction mechanism of lo- 
production. The orientation of the oj-spin is correlated 
with that of the decay plane of the u! meson. The lat- 
ter can be determined from the pion momenta in the 
CMS of the decaying oj-meson. In our experiment the 
two protons and the charged pions are detected, however 
the momentum resolution for these fast pions (^ « 1) 
is not sufficient to make direct use of the measured mo- 
menta such as to calculate the missing 4-momentum of 
the tt". Nevertheless, its 4-momentum can be determined 
by exploiting the time of fiight of the two observed pions, 
TO-Fobs,Tr, and the w- momentum. The latter is known as 
Puj — ?'(pPinitiai)~-P(PPfinai) ^ud the 4-momenta of all three 
pions must add up to Puj = Ptt- + P-k+ +P-kO- By varying 
the 4-momenta of the charged pions and calculating the 
corresponding time of fiight, T0i^caic,7r, one can find that 
momentum partition which yields the minimum deviation 
{ATOFf = {TOF,,i,^^, - TOF,i,,,^,y + {TOF,,i,,^, - 
TOFobs.TTa)^- This method allows to deduce all pion mo- 
menta and thus the orientation of the w decay plane which 
is determined with an accuracy of cr = 10°. We like to 
point out, however, that this analysis cannot be used to 
improve on the w-signal as the background below the sig- 
nal cannot be reduced. 

It was pointed out by Gottfried and Jackson [33] and 
later on by Titov et al. [TB], [3S] that it is essential to 
speciiy the reference frame within which an observable 
is measured. These authors showed that for the reaction 
PP PP'/'j 4> ~^ K+K~ the decay angle 9, defined as the 
polar angle of the direction of flight of one of the decay 
particles in the 0-meson's rest frame with respect to the 



beam direction, shows an angular distribution (normalised 
to 1) given by 

WicosO) ^ ^(1 - Poo + (3poo - l)cos20) (3) 

If 0-production is considered just above threshold the 
spin of the 0-meson is aligned and the spin density matrix 
element poo = 0; then eq.Q reduces to 

W{cose) = ^sin^e (4) 

This relation was used by e.g. Balestra et al. [3 , Rekalo 
et al. and Hartmann et al. (37) in their investigations 
of the spin alignment of the 0-meson. 

The appropriate reference axis for the 3-body decay of 
the w-meson is the normal to the decay plane. The angle 
of the normal with respect to the beam axis becomes 7 = 
tt/2~9 and replaces the angle 9 in eg.® when translating 
the 2-body 0-decay into the 3-body w-decay. Eq.(|3|) then 
reads 

W^(cos7) ^ N(l - Poo + (3/500 - l)sm^l) 

= N(2poo - (3poo - l)cos27) (5) 

with the proper normalisation constant N = l/(2poo+2/3). 

We will exploit the eq.® in order to investigate a pos- 
sible spin alignment in the w-production. The term align- 
ment arises from the following consideration: the orbital 
angular momentum of the initial state of the pp — ^ ppw re- 
action is always oriented perpendicularly to the beam and 
the projection quantum number of the total angular mo- 
mentum Mj. of the pp-system can only assume the values 
of or ±1 as the Clebsch-Gordan coefficients show. The 
restriction Mj. = 0, ±1 carries over to the final state, and 
the notion of an aligned total angular momentum applies 
if Mj^ = ±1 only. 

Directly at threshold, the transition ^Pi ^SqS (am- 
plitude /i of table [2]) is the only possible one. In this case, 
Ji = Jf = Sui = i and Mj. = nis^ = ±1, i-e. the spin of 
the cj-meson is aligned (which goes along with poo = 0) 
and the angular distribution becomes according to eq. ([5]) 

3 

W{cosj) = -cos^7 (6) 

This finding, solely based on conservation laws which 
are reflected by the appropriate Clebsch-Gordan coeffi- 
cients, is true independent of the reaction model under 
consideration. 

It should be pointed out, that poo holds only at 
threshold but assumes values of about 0.2 already slightly 
above. This is a consequence of the spin-orbit interaction 
[38] . Naturally, this effect increases with the vector meson 
momentum q. 

Above threshold, we have to consider flrstly, Sp type 
final states (amplitudes /2 and /a of table [3]). For the am- 
plitude /2 the initial state ^Sq does not exhibit any orien- 
tation, hence the w-spin has no particular orientation and 
W{cosj) = 1/2 is isotropic. This reflects the fact that all 
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Table 5. Parameters of the function eq.lO used to describe the 
distribution of the angle 7 between the normal on the oj-decay 
plane and the beam-axis. 



£ (MeV) 




poo 


92 


8.3 ± 1.2 


0.17 ±0.07 


128 


12.0 ± 0.6 


0.19 ±0.03 


173 


28.0 ±2.3 


0.24 ±0.05 



spin projection states are equally populated and, hence, 
the spin density matrix element poo — Pii = P-i-i = l/3. 
For the amplitude /a (^D2 — > ^Sop) one finds Mj^ = 0, 
i.e. the total angular momentum of the w- meson is perpen- 
dicular to the beam axis. The w-spin projection quantum 
number, however, can assume the values rris^ = 0,±1 re- 
sulting in a distribution of the type a ■ sin^7 + b ■ cos^j 
which tends to reduce the anisotropy due to amplitude fi . 
Secondly, we consider Sd type final states. Distributions of 
the same type as for amplitude /a are expected. Thirdly, if 
we take into account Ps type final states we find isotropy 
for amplitude f4, type a ■ sin^7 ± b ■ cos^j distributions for 
fs and alignment for fg. 

The measured angular distribution of the w-spin direc- 
tion is shown in fig. [S] We describe these angular distri- 
butions by employing eq.([5]) 



da 



dcos7 2/900 ±2/3 



(2poo- (3poo- l)cos27). (7) 



Fig. [5] depicts the results of the fitting as a solid line. 
The values determined for atot and (poo) are listed in ta- 
ble [3 

The increasing value of poo indicates the transition 
from an aligned w-spin (poo = 0) to an arbitrarily oriented 
w-spin (poo — 1/3). This is due to two effects, firstly the 
increasing importance of reaction amplitudes that show 
no alignment, namely those of the initial state "'^Sq which 
provide an isotropic contribution to the angular distribu- 
tion as well as those which feature a sin^7-dependence, 
and secondly the the above mentioned increase of poo due 
to spin orbit interaction being important above threshold 
[55] . The energy dependence of the parameter poo sets a 
benchmark for theoretical models of w-production. 



3.4 Invariant mass spectra 

Invariant mass spectra of the two-body subsystems can 
be used to search for deviations from phase space which 
would be indicative of a resonance. In order to determine 
an invariant mass spectrum the w-signal needs to be ex- 
tracted from the missing mass spectra generated for vari- 
ous bins in the invariant mass. It turned out that this was 
possible only if a Monte Carlo simulated pion-cocktail was 
used for the determination of the background. The statis- 
tical basis of our data does not allow to generate a Dalitz 
plot. 



The results obtained are shown in figs. [TU] and [TT] (up- 
per part of each frame, the lower part shows the accep- 
tance) together with the simulated distribution based on 
an equally populated phase space for the w-production. 
The fuzzy fringes of the simulated distributions (upper 
end of Mpp and lower end of M^p) reflect the finite width 
of the w-meson; the experimental data are reproduced very 
well. A^* -resonances with widths in the order of 100 to 200 
MeV do not cause the invariant mass spectrum to deviate 
perceptibly from phase space, as simulations have shown. 
Arguing the other way around we conclude that any res- 
onance with a width below our missing mass resolution 
of 10 MeV would show up if their cross section would be 
above w 0.5pb. 

It is obvious from these results that indications of reso- 
nant w-production are missing. This finding is compatible 
with that from the investigation of the proton angular dis- 
tributions. An observable closely related to the invariant 
mass will be discussed in the following subsection. 

3.5 Helicity and Jackson angle 

3.5.1 Notation 

In a 2 — > 3 reaction {a + b ^ 1±2±3) with conservation 
of 4-momentum given by 



Pa + Pb ^ Pi + P2 + P3 



(8) 



the initial state in the overall CMS is defined by a direction 
and the total energy ^/s which both arc known in standard 
collision experiments. The description of the final state 
requires 3x3 = 9 variables (besides the masses in the 
exit channel) which, due to conservation of momentum 
and energy, reduce to five variables. Without polarisation 
in the entrance channel the azimuthal orientation of the 
final state is arbitrary, reducing the number of necessary 
variables to four. 

The four variables needed to uniquely describe the re- 
action can be represented by two invariants and two an- 
gles. For the latter one may choose e.g. the CMS angles of 
any two final state particles; for the invariants one often 
chooses the invariant masses of 2-body subsystems 

(9) 
(10) 



Sl3 = [Pl + Pj.f = {Pa +Pb~ P2f 
S23 = {P2 ± Pzf {Pa +Pb- Plf 



where S12 may be chosen as an alternative to S13 or S23- 
With this choice, only quantities of the exit channel are 
taken into account. This shortcoming can be avoided by 
choosing a suitable angle in an appropriate Lorentz refer- 
ence system [HI]- The latter is obtained by boosting the 
final state in such a way that the momenta of particle 2 
and 3 add up to zero; one then obtains from eq. (jH]) 

= P2 ± P3 = Pa ± {Pb - Pl) =Pb+ {Pa - Pl)- (f 1) 

Now, an angle connecting exit and entrance channel 
is the (polar) Jackson angle (notation of Byckling and 
Kajantie [39j) 



cosC' 



Pb P3 



\Pb\ ■ \P3\ 



P2 = 



-P3 



(12) 
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Fig. 9. Distributions of the angle between the direction of the iLj-spin and the beam axis at the three excess energies (upper 
part in each frame; lower part: acceptance). See text for details. 
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Fig. 10. Invariant pp-mass distribution (upper part of each frame, acceptance: lower part). The histograms represent the phase 
space distributions. 
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Fig. 11. Invariant tjp-mass distribution (upper part of each frame, acceptance: lower part). The histograms represent the phase 
space distributions. 



The superscript specifies the Lorentz reference system R23 
chosen through p2 = —ps', the first subscript indicates 
which particle defines via its momentum direction the po- 
lar reference axis (a or b for the Jackson frame). The sec- 
ond subscript specifies which particle is used when calcu- 
lating the Jackson angle. Since p2 = — Ps, the choice of 
this particle is arbitrary, and 9f^^ — tt ~ 0^^^. 

In the same Lorentz reference system one can also cal- 
culate the (polar) helicity angle 

cosefi' = ,^Vf I • (13) 



superscript, while the first subscript indicates which par- 
ticle defines via its momentum direction the polar axis 
(1 for the helicity frame). The second subscript specifies 
which particle is used when calculating the helicity angle. 
Since P2 — — Pa, the choice of this particle is arbitrary, 
and 9f-^^ — tt — 0^2^. Of course, all formulae are also valid 
for any other Lorentz reference systems (superscript R12 
or R13 in eq. p^ or ea.p^ with the corresponding first 
and second subscript. 



3.5.2 Helicity angle 



which, in contrast to the Jackson angle depends on exit 
channel properties only. Again the Lorentz reference sys- 
tem "R23" chosen through p2 = — ps is specified by the 



The helicity angle relates particles 2 and 3 to particle 1, 
and hence carries information solely about the exit chan- 
nel. This information is also contained in the Dalitz-plot: 
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Fig. 12. Helicity angle distribution of the w-meson in the pp-helicity frame (upper part of each frame, acceptance: lower part). 
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Fig. 13. Helicity angle distribution of the "other" proton in the wp-helicity frame (upper part of each frame, acceptance: lower 
part). 
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Fig. 14. Distribution of the angle between the beam direction and that of a proton in the pp-Jackson frame (Jackson angle) 
(upper part of each frame, acceptance: lower part). 
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Fig. 15. Distribution of the angle between beam and the w-meson in the wp- Jackson frame (Jackson angle) (upper part of each 
frame, acceptance: lower part). 
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A resonance in the 2-3-systeni is observed as an enhance- 
ment of events along the S23-axis at a certain mass and 
the distribution of cos 9^^^ is simultaneously seen along 
the Si3-axis. From this helicity angle distribution the spin 
of the resonance can be inferred in particular cases. It 
should be noted that this angular distribution is not nec- 
essarily symmetric with respect to cos 9^^^ = as it does 
not relate to the entrance channel. 

Two helicity angles can be considered for the reaction 
pp ppt^: In the first case the two protons are used to 
define the Lorentz reference system and the polar helicity 
angle of the w-meson is determined with respect to one of 
the protons. Since the two protons are indistinguishable 
we must take into account two distributions, namely 9^pp 
and6'J}PP. This results in a symmetrised helicity angle dis- 
tribution. 

The result is shown in fig. [12] where all possible values 
of Si2 are accounted for. This corresponds to a summation 
of Dalitz-plot entries with constant values of cos 9^^^. The 
solid lines represent the results of least square fitting using 
eq.©. The coefficients are listed in the first section of 
table El 

At e = 92 MeV one observes isotropy as well as at e = 
128 MeV. At £ = 173 MeV the increase of cross section 
at cos^J'p^ = ±1 seems to indicate that the w-meson has 
some tendency to move along with a proton. This may 
reflect the angular momentum of the w-meson with respect 
to the ^Sq di-proton and not necessarily a contribution of 
a resonance. A theoretical model is asked for to explain 
this observation. 

In the second case one proton (for instance particle 
2) and the w-meson (particle 3) are used to define the 
Lorentz reference system, and the polar helicity angle of 
the other proton (particle 1) is determined with respect 
to the direction of the w-meson. Since the two protons 
are indistinguishable we must average the distributions of 
(^p^^S^ 9,nd 9^^J'^ for a given event. A possible anisotropy 
in the helicity angle is not destroyed by this procedure. 

Fig. [13] shows this helicity angle distribution where all 
possible values of S23 were taken into account. This corre- 
sponds to a summation of Dalitz-plot entries with constant 
values of cos9^p. The solid lines represent the results of 
least square fitting using eq.Q but allowing Pi in addition 
with a weight of a'l — ai/ao. The coefficients are listed in 
the second section of table [5] 

The helicity angle distributions at e = 92 and 173 MeV 
are isotropic within uncertainties. At e = 128 MeV sig- 
nificant anisotropy and asymmetry values are observed 
pointing at an inhomogeneous population of the Dalitz 
plot which is probably caused by the angular momentum 
in the exit channel. This finding, however, is not in contra- 
diction with the helicity angle distribution of fig. as a 
different projection of the Dalitz plot is presented. Again, 
without a theoretical model one cannot explain the origin 
of this effect. 



3.5.3 Jackson angle 

The Jackson angle relates the direction of the beam mo- 
mentum with the orientation of the axis given by p2 = 
— P3. Ea.([TT]) can be interpreted, as suggested by Gottfried 
and Jackson [33], as a 2 — >■ 2 reaction, namely 2-1-3 — > b+x 
where particle x with Px = Pa — Pi can be identified with 
the exchange meson in diagram (a) or (c) in fig.[T] Thus it 
is clear that any structure in the Jackson angle distribu- 
tion gives direct information on the angular momentum of 
the system R23 which could, but not necessarily has to, 
be a resonance [3D]. 

Two Jackson angles can be considered for the reaction 
pp — > ppw: In the first case the two protons (particle 1 
and 2) are used to define the Lorentz reference system 
and the angle between the direction of one of the two pro- 
tons with respect to the beam direction is called the polar 
Jackson angle. Since the protons are indistinguishable we 
must take into account the angles 9^^^ and ^^p^'', a proce- 
dure which leads to a distribution of a symmetrised polar 
Jackson angle 9^^^. Since beam and target proton can 
not been distinguished, both Jackson angles with respect 
to beam and target are taken into account for each event 
(maintaining for simplicity the sudscript b in 9^^^ ) 

In fig. [33] we show this symmetrised Jackson angle dis- 
tribution. The solid lines represent the results of least 
square fitting using eq.©. The coefficients are listed in 
the third section of table [H] 

The fits suggest a slight anisotropy which implies that 
there is at most a very weak correlation between the beam 
direction and the final proton pair. This is compatible 
with the proton angular distributions presented in fig. [3] 
which showed an insignificant contribution of Ps type final 
states. 

In the second case one proton (for instance particle 1) 
and the w-meson (particle 3) are used to define the Lorentz 
reference system, and the polar Jackson angle is measured 
as the direction of the w-meson with respect to the beam 
direction. Since the two protons in the exit channel are 
indistinguishable we must average the angles 9^^^ and 
^huj^^ in order to yield 9^^. This averaging causes a di- 
lution of a signal from a potential resonance since this is 
either found in the wpi or ajp2 system, while the non- 
resonant one furnishes an uncorrelated, hence isotropic 
background. 

Again, the indistinguishability of beam and target are 
taken into account. However, in the Lorentz reference sys- 
tem beam and target momenta are not aligned, hence this 
procedure is not a symmetrisation, but the distribution 
turns out to be symmetric. 

In fig.[Tn]we show this Jackson angle distribution. The 
solid lines represent the results of least square fitting using 
eq.((2]). The coefficients are listed in the fourth section of 
table [S] 

The Jackson angle distribution at e = 92 MeV is al- 
most isotropic but pronounced anisotropics are observed 
at the other excess energies beyond doubt. This shows that 
there is a non-zero relative angular momentum in the up 
system. However, a conclusion that this is a direct signal 
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Table 6. Coefficients of Legendre polynomials determined by fitting eq.((2l) to the distributions of the different observables. 



observable 


£ (MeV) 


(Ttot (/ib) 




a'2 


cos e^i;^ 

(fig.[l2j 


92 
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173 


8.8 ± 0.5 

12.7 ±0.3 

26.8 ± 1.2 




0.01 ± 0.14 
-0.08 ±0.07 
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COS f7p^ 

(fig. [SI 


no 
yz 

128 

173 


y.i ± U.O 

12.7 ±0.3 

26.8 ± 1.2 


n A? 1 nil 

— U.Ul ± U.ii 

0.07 ±0.04 
-0.04 ±0.08 


n AQ 1 Air 

U.Uo ± U.iO 

-0.17 ±0.06 
-0.06 ±0.11 


cos<7 

(fig. [a 


92 
128 
173 


9.2 ±0.5 
12.7 ±0.2 
27.5 ±0.9 




-0.13 ±0.15 
0.09 ± 0.05 
0.22 ± 0.09 


cose^^P 

(fig. [13 


92 
128 
173 


8.7 ±0.4 
12.7 ±0.2 
26.5 ±0.9 


-0.01 ±0.08 
-0.04 ± 0.04 
0.01 ±0.06 


0.12 ±0.11 
0.35 ± 0.04 
0.37 ±0.08 




0.5 „ 1 

cose7 



0.5 „ 1 
cose^r 



Fig. 16. Comparison of the experimental results with a Monte Carlo simulation using the measured distribution of cos6'„ as a 
weight function. 



of a iV*-resonance decaying into pcj should be drawn with 
caution. The deduced values of L — 0,1,2, as discussed 
in the context of the uj angular distributions (fig. [7]), are 
reflected by the Jackson angular distribution. In order to 
corroborate this argument we present in fig. (together 
with data already shown) results of Monte Carlo simu- 
lations at e = 128 MeV where the event generator was 
modified such as to reproduce the w angular distribution 
of fig. [71 again shown in the upper left frame of fig. 1161 
This weight function on the cj angular distribution also 
modifies the proton angular (lower left) and Jackson an- 
gle distributions (upper middle and right frame) in such 
a way that they perfectly match the experimental ones. 
We obviously observe in the Jackson angle distribution 
the deviation of the reaction kinematics from pure phase 



space. The helicity angle distributions (lower middle and 
lower right frame), however, cannot be reproduced. 

We tentatively conclude that we observe the influence 
of some particular reaction dynamics: the transition ma- 
trix element connecting the initial with the final state 
shows a dependence on q and but not on p or L, creat- 
ing anisotropic angular distributions of the w-mesons, he- 
licity angle 9^^^, and Jackson angle O^'^. The final state 
protons are in the ^Sq state, resulting in isotropic angular 
distributions of 0*, helicity angle 0^^^, and Jackson angle 

We would like to note that the results of the corre- 
sponding investigation at e = 173 MeV allow the same 
conclusion. 
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4 Summary 

In this paper we presented a systematic study of the pro- 
duction of cj-mcsons in proton-proton-collisions, carried 
out in a kincmatically complete experiment at three excess 
energies of £ = 92, 128, 173 MeV. The large-acceptance 
COSY-TOF spectrometer allows the almost unambigu- 
ous and simultaneous identification of different reaction 
channels. We described in detail how the yield of w-evcnts 
was determined in the presence of an unavoidable physical 
background which is due to the production and decay of 
the p-meson [p — > tt+tt^) and non-resonant tt+tt" as well 
as Tr+TT^Tr" production. Total cross sections, angular dis- 
tributions both of w-mesons and protons were measured 
in the overall CMS as well as helicity and Jackson angle 
distributions in both the pp and wp helicity and Jackson 
frames, respectively. In addition, the orientation of the lo- 
spin and invariant mass spectra were determined. All total 
cross sections obtained by integrating the various differ- 
ential distributions agree within uncertainty, proving the 
consistency of our data evaluation. 

The major results are as follows: We conclude from 
the angular distribution of the w-meson that its produc- 
tion takes place dominantly in Ss and Sp final states for 
the lower energies, where only little room is left for Ps final 
states, and, additionally, type Sd final states for the high- 
est energy. The production of w-mesons via A'^*-resonances 
was ruled out to be the major reaction mechanism. It is, 
however, conceivable that resonant w-production via the 
broad subthreshold resonances Sn and D13 as well as res- 
onances above threshold, Pn and P13, may happen for 
particular initial states. However it probably plays a mi- 
nor role among the various lo production mechanisms dis- 
cussed in the literature. Invariant mass spectra for both 
the pp and pw subsystems are found to be compatible 
with phase space distributions underlining the minor im- 
portance of resonant w-production. A dominant role of 
^Pi and ^So initial partial waves for w-production was 
concluded from the orientation of the w-spin. Although 
we observe anisotropic Jackson angle distributions in the 
wp- Jackson frame we argue that this is not an indication 
of a resonance but rather a kinematical effect reflecting 
the anisotropy of the w angular distribution in the over- 
all CMS. The helicity angle distribution in the wp- helicity 
frame shows an anisotropy which, in addition to the ori- 
entation of the w-spin, is probably the most sensitive ob- 
servable to judge the validity of various theoretical de- 
scriptions of the production process. 
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Table 7. Differential cross sections in yub/sr of the angular distributions. The central value of a cosine interval of 0.2 or 0.1 is 
listed. 
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0.05 


1.80 ±0.53 


2.16 ±0.35 


1.61 ± 0.46 


2.12 ± 0.37 


2.13 ±0.39 


1.30 ±0.24 




0.15 


2.37 ± 0.49 


2.42 ± 0.49 


2.05 ± 0.41 


2.21 ±0.37 


1.83 ± 0.32 


1.81 ±0.31 




0.25 


2.28 ± 0.39 


2.23 ± 0.38 


2.10 ±0.45 


2.31 ±0.39 


1.56 ± 0.36 


1.64 ±0.31 




0.35 


1.96 ±0.31 


2.90 ±0.51 


2.30 ± 0.43 


2.85 ±0.41 


2.17 ±0.23 


1.78 ± 0.24 




0.45 


2.41 ±0.40 


1.78 ±0.51 


2.04 ± 0.39 


2.34 ± 0.48 


2.25 ±0.32 


2.08 ± 0.40 




0.55 


1.87 ±0.29 


1.92 ±0.50 


2.25 ±0.51 


1.77 ±0.44 


2.21 ±0.32 


1.89 ± 0.29 




0.65 


2.24 ± 0.35 


2.00 ± 0.66 


2.13 ±0.44 


1.82 ±0.43 


2.07 ± 0.29 


2.37 ±0.30 




0.75 


2.27 ± 0.39 


3.51 ± 0.81 


1.75 ± 0.38 


2.14 ±0.40 


2.47 ± 0.33 


2.00 ± 0.38 




0.85 


3.04 ± 0.36 


2.47 ± 0.78 


2.28 ± 0.47 


2.27 ± 0.43 


2.51 ± 0.33 


2.37 ±0.38 




0.95 


4.04 ± 0.41 


2.06 ± 1.24 


2.86 ± 0.37 


1.49 ± 0.43 


2.41 ± 0.57 


3.29 ± 0.33 
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Table 8. DiflFerential cross sections in fib/sT for the orientation 
of the uj decay plane at the different excess energies. The central 
value of a cosine interval of 0.2 is listed. 



cos 7 


e = 92 MeV 


e = 128 MeV 


s = 173 MeV 


0.1 


0.58 ± 0.45 


0.65 ± 0.09 


2.21 ± 0.57 


0.3 


0.35 ±0.17 


0.76 ± 0.06 


1.78 ± 0.48 


0.5 


0.87 ±0.18 


0.96 ± 0.04 


2.13 ± 0.28 


0.7 


0.70 ±0.15 


1.04 ±0.03 


2.21 ±0.45 


0.9 


0.94 ±0.15 


1.28 ±0.03 


2.85 ± 0.33 



Table 9. Differential cross sections in /ib/sr for the invariant 
mass of the wp-system, statistical errors are given only. The 
central mass of an interval of 10 MeV/c^ is listed. 



M„p (MeV/c^) e = 92 MeV e = 128 MeV e = 173 MeV 



1715 


0.16 


± 


0.07 


0, 


.16 


± 


0. 


.01 


0. 


.30 


± 


0. 


.07 


1725 


0.45 


± 


0.08 


0, 


.40 


± 


0. 


.01 


0. 


.54 


± 


0. 


.10 


1735 


0.91 


± 


0.11 


0, 


.78 


± 


0. 


.02 


1. 


.21 


± 


0. 


.17 


1745 


1.18 


± 


0.12 


1, 


.20 


± 


0. 


.02 


1. 


.58 


± 


0. 


.24 


1755 


1.28 


± 


0.12 


1. 


.19 


± 


0. 


.02 


1. 


.34 


± 


0. 


.22 


1765 


1.17 


± 


0.13 


1, 


.31 


± 


0. 


.02 


1. 


.97 


± 


0. 


.24 


1775 


1.26 


± 


0.13 


1, 


.27 


± 


0. 


.02 


1. 


.65 


± 


0. 


.24 


1785 


1.37 


± 


0.13 


1, 


.27 


± 


0. 


.02 


2. 


.52 


± 


0. 


.25 


1795 


1.24 


± 


0.11 


1, 


.28 


± 


0. 


.02 


2. 


.07 


± 


0. 


.24 


1805 


0.78 


± 


0.10 


1, 


.13 


± 


0. 


.02 


1. 


.89 


± 


0. 


.23 


1815 








1, 


.25 


± 


0. 


.02 


2. 


.08 


± 


0. 


.23 


1825 








0, 


.99 


± 


0. 


.02 


1. 


.67 


± 


0. 


.23 


1835 








0, 


.64 


± 


0. 


.02 


2. 


.26 


± 


0. 


.23 


1845 








0. 


.42 


± 


0. 


.02 


2. 


.11 


± 


0. 


.22 


1855 


















1. 


.74 


± 


0. 


.20 


1865 


















1. 


.81 


± 


0. 


.20 


1875 


















1. 


.49 


± 


0. 


.18 


1885 


















0. 


.86 


± 


0. 


.15 



Table 10. Differential cross sections in /ib/sr for the invariant 
mass of the pp-system, statistical errors are given only. The 
central mass of an interval of 10 MeV/c^ is listed. 



Mpp (MeV/c^) 


e = 92 MeV 


e = 128 MeV 


e = 173 MeV 


1882 


0.95 ± 0.15 


0.69 ± 0.02 


1.27 ±0.23 


1892 


1.05 ±0.15 


0.77 ± 0.02 


1.29 ±0.24 


1902 


1.29 ±0.17 


1.10 ±0.03 


0.83 ±0.26 


1912 


1.20 ±0.18 


1.03 ± 0.03 


1.07 ±0.25 


1922 


1.32 ±0.18 


1.24 ±0.03 


1.73 ±0.28 


1932 


1.44 ±0.18 


1.59 ± 0.03 


2.77 ± 0.29 


1942 


0.94 ±0.15 


1.41 ± 0.03 


1.99 ± 0.35 


1952 


1.25 ±0.17 


1.34 ±0.03 


2.40 ± 0.31 


1962 


0.41 ±0.15 


1.30 ± 0.03 


2.48 ±0.33 


1972 


0.01 ±0.16 


1.11 ±0.03 


1.70 ± 0.32 


1982 




0.97 ±0.03 


2.04 ± 0.31 


1992 




0.44 ± 0.02 


2.60 ± 0.31 


2002 




0.23 ±0.01 


2.07 ±0.32 


2012 






2.00 ±0.33 


2022 






1.55 ±0.30 


2032 






0.89 ±0.26 


2042 






0.91 ±0.25 



